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MIXED  MODE  FRACTURE  OF  CONCRETE 
AT  HIGH  STRAIN  RATE 

Mixed-mode  fracture  specimens  were  tested  at  quasi-static  and  impact  rate  of 
loading.  Laser  holographic  interferometry,  reflection  photoelasticity,  image  processing 
and  high  speed  photography  were  used  for  recording,  measurement  and  analysis  of  data. 

Center  notched  plate  specimens  were  loaded  in  tension.  Two  dimensional 
displacement  fields  were  measured  during  stable  crack  propagation  using  holographic 
interferometry.  Image  processing  was  used  to  evaluate  the  fringe  patterns.  The  existence 
of  forces  transmitted  through  the  crack  faces  was  associated  with  tensile  strain  behind 
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fracture  was  used  to  predict  experimental  observations. 

Holographic  interferometry  was  generalized  to  measure  three  dimensional 
displacement  fields.  Strain  fields  and  crack  growth  were  measured  for  a  diametral  disk 
specimen  (Brazilian  Test).  Tensile  strain  levels  before  cracking  in  the  Brazilian 
configuration  are  significantly  higher  than  those  observed  in  the  center  notched  plates. 
Crack  initiation  occurs  at  approximately  70%  of  the  ultimate  load.  Cracks  do  not  initiate 
in  the  region  of  uniform  tensile  stress  as  is  commonly  assumed  in  calculating  tensile 
strength.  Collapse  of  the  specimen  is  due  to  the  reported  ’wedge’  formation  and  occurs 
after  a  single  crack  has  propagated  through  the  specimen  and  a  significant  loss  of 
strength  has  taken  place. 

Center  notched  disk  specimens  were  tested  under  diametral  loading.  By  rotating 
the  notch  angle  with  respect  to  the  loading  direction,  it  was  possible  to  change  the  failure 
mode  from  mode  I  to  mixed  mode.  The  size  of  the  FPZ  defined  as  a  function  of  the 


tensile  strain  distribution  depends  on  the  fracture  mode.  The  apparent  stress  intensity 
factor  showed  a  R-curve  type  of  behavior. 

Holographic  observations  were  correlated  with  measurements  obtained  using 
reflection  photoelasticity  at  quasi-static  loading  for  the  center  notched  disk  specimen.  To 
study  the  influence  of  rate  of  loading,  an  instrumented  Charpy  type  of  test  system  was 
modified.  Center  notched  disk  specimens  subjected  to  impact  loading  were  studied  with 
reflection  photoelasticity  and  high  speed  photography  (10,000  fps).  The  increase  in 
fracture  toughness  at  crack  initiation  for  high  rates  of  loading  is  affected  by  the  failure 
mode.  The  increase  in  apparent  fracture  toughness  for  tensile  cracking  is  approximately 
double  the  corresponding  value  for  shear  cracking.  Failure  of  both  quasi-static  and 
impact  rate  of  loading  was  characterized  by  crack  initiation  and  stable  crack  growth.  After 
initial  crack  growth,  the  crack  continues  to  propagate  in  mode  I  regardless  of  the  notch 
angle. 

Reflection  photoelasticity  can  be  used  to  determine  crack  initiation  and  crack 
growth  at  both  quasi-static  and  impact  rate  of  loading.  However  it  is  not  sensitive  enough 
to  provide  information  regarding  the  state  of  deformation  away  from  the  crack  tip. 
Holographic  interferometry  with  multiple  sensitivity  vectors  and  digital  image  analysis  can 
provide  a  three  dimensional  displacement  field  with  an  accuracy  of  a  1/4  micron. 

To  design  reinforced  concrete  structures  subjected  to  impact  loading,  one  must 
consider  the  observation  that  the  influence  of  rate  of  loading  is  different  for  flexural  failure 
and  shear  failure. 
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SECTION  I 


INTRODUCTION 


A.  OBJECTIVES 

The  purpose  of  the  investigations  reported  in  this  document  is  to  obtain  fracture 
parameters  for  concrete  subjected  to  mixed  mode  type  loading  at  varying  strain  rates 
(ICT6  to  1  per  sec).  Center-notched  disk  specimens  (CNDS)  were  tested  at  slow  rates  in 
a  closed-loop,  servo-controlled  test  machine  and  at  intermediate  rate  of  loading  (strain 
rate  1  to  10  per  sec)  in  an  instrumented  Charpy  machine.  This  investigation  will  be 
complemented  by  high  rates  of  loading  (strain  rate  >  103  per  sec)  obtained  using  a  SHPB 
setup  at  the  Air  Force  Engineering  Service  Center  at  the  Tyndall  Air  Force  Base  (Figure 
Figure  1).  The  studies  at  slow  rates  were  conducted  to  investigate  the  development  of 
fracture  process  zone  in  Mode  I  and  mixed  mode  using  laser  holographic  interferometry. 
Additionally,  laser  holographic  interferometry  was  used  to  calibrate  reflection 
photoelasticity  recordings  on  similar  specimens  at  comparable  rate  of  loading.  The 
experimental  techniques  developed  at  slow  rates  were  then  extended  to  impact  rates  of 
loading.  For  impact  loading,  crack  paths  were  traced  by  reflection  photoelasticity. 

B.  BACKGROUND 

Experimental  results  indicate  that  the  tensile  and  compressive  strengths  of  concrete 
increase  with  increase  in  rate  of  loading  (References  1  and  2).  There  is  an  observed 
dependence  between  the  failure  mechanism  and  the  rate  of  loading  (Reference  3).  A 
structure  may  fail  in  a  ductile  bending  mode  at  slow  rates  of  loadings  and  in  shear  at 
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higher  rates  of  loading.  To  predict  the  rate  dependent  failure  mechanisms  in  reinforced 
concrete  structures,  the  knowledge  of  the  effect  of  impact  loading  on  fracture  of  concrete 
is  essential.  In  the  most  general  case,  the  fracture  of  concrete  occurs  under  combined 
fracture  modes.  However,  little  information  is  available  on  the  effects  of  rates  of  loading 
on  fracture  properties  of  concrete  subjected  to  combined  tensile  and  shear  loading. 

Recently,  strain  rate  effects  on  compressive  and  tensile  strength  of  concrete  were 
conducted  by  Ross,  et.al.,  Malvern,  et.  al.  and  Reinhardt,  et.  al.  (References  4,  5  and 
6)  using  Split  Hopkinson  Pressure  Bar  (SHPB)  tests.  Ross  (Reference  4)  used  concrete 
cylinders  loaded  in  uniaxial  compression  and  in  diametral  compression  (Brazilian  Test)  to 
study  compression  and  tension  at  high  strain  rates.  By  introducing  a  central  notch  on  the 
disk  specimen,  mixed  mode  fracture  conditions  can  be  obtained  with  a  wide  range  of 
proportions  between  the  tensile  and  shear  components  (Figure  2). 


Figure  2  K,/KIC  vs.  Kn/Kllc  Curve  for  Different  Load  Configurations 


analysis  and  how  they  can  be  used  to  evaluate  two-dimensional  strain  fields  and  how  'full 


filed  displacement  measurements  can  be  used  to  characterize  the  fracture  process.  In 
Section  III,  an  extrapolation  of  the  techniques  described  in  Section  II  for  three  dimensional 
displacement  field  evaluation  under  closed  loop  loading  conditions  are  described.  Using 
this  techniques,  crack  initiation  and  propagation  was  evaluated  in  diametrically  loaded 
cylinder  specimens.  Section  IV  presents  the  results  of  the  strain  field  evaluation  on  CMOS 
at  different  notch  inclination  angles.  Experimental  observations  were  associated  with  the 
mode  of  fracture,  in  Section  V  photoelasticity  recordings  at  slow  rate  of  loading  are 
compared  with  the  holographic  recordings  at  slow  rates  of  loadings.  An  experimental 
arrangement  for  high-speed  photoelasticity  to  be  used  for  impact  loading  is  presented 
and  the  results  are  compared  with  the  slow  rate  of  loading  evaluations. 


SECTION  II 


STRAIN  FIELDS  IN  MODE  I  FRACTURE  PROCESS  ZONE 


A.  INTRODUCTION 

The  fracture  of  concrete  is  characterized  by  formation  and  propagation  of  fracture 
process  zone.  Several  phenomena  have  been  associated  with  fracture  process  zone 
including  microcracking  around  the  crack  tip,  localization  and  strain  softening  and  crack 
bridging  (often  also  termed  aggregate  interlock).  A  key  to  a  better  understanding  of 
fracture  mechanics  of  concrete  is  the  accurate  observation  of  fracture  process  zone  on 
a  microscopic  scale.  The  research  in  this  study  has  been  directed  to  better 
understanding  of  the  role  of  fracture  process  zone. 

The  first  observation  of  this  kind  was  made  by  Cedolin,  Dei  Poli  and  lori  (Reference 
7).  They  used  laser  Moire  interferometry  for  wide  field  measurements  in  a  specimen 
under  tension.  Such  measurements  allowed  an  indirect  determination  of  the  local  stress- 
strain  and  stress-crack  opening  displacement  (COD)  relationships.  White  light  Moire 
interferometry  was  also  used  by  Du,  Kobayashi  and  Hawkins  (Reference  8)  to  measure 
crack  opening  displacement  along  the  fracture  process  zone  (FPZ)  for  concrete 
specimens  under  three  point  bending.  Regnault  and  Bruhwiler  (Reference  9)  used  Moire 
holographic  interferometry  to  measure  strain  field  on  wedge-splitting  test  specimens  and 
concluded  that  the  model  based  on  the  exact  knowledge  of  traction  free  crack  tip  location 
does  not  describe  the  material  behavior  in  a  realistic  way.  Raiss,  Dougill  and  Newman 
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(Reference  10)  also  used  laser  Moire  interferometry  to  measure  strain  field  on  direct- 
tension  specimens. 

The  accuracy  of  displacement  observation  with  laser  Moire  interferometry  is  limited 
by  spacing  of  grating  and  possible  debonding  of  grating.  Miller,  Shah  and  Bjelkhagen 
(Reference  11)  used  laser  sandwich  hologram  interferometry  to  obtain  crack  opening 
displacements  with  an  accuracy  of  1/4  of  a  micrometer.  With  their  technique  it  was  not 
possible  to  measure  wide-field  two-dimensional  displacement  fields  (only  crack  opening 
displacements  were  able  to  be  measured)  since  the  loaded  specimen  was  illuminated  by 
a  laser  beam  from  only  one  direction.  In  the  technique  described  in  this  report,  three 
simultaneous  holograms  were  taken  by  illuminating  the  specimen  from  three  different 
directions.  In  addition,  digital  image  analysis  techniques  were  developed  to  accurately 
evaluate  wide  field  in  plane  deformation.  The  results  of  the  strain  field  calculation 
obtained  using  this  newly  developed  technique  for  mortar  center- notched  plate  specimens 
are  discussed  in  this  investigation. 


B.  EXPERIMENTAL  PROGRAM 

Center-notched  mortar  plate  specimens  were  used  (Figure  3).  One-inch  (25  mm) 
diameter  holes  with  3/8  inches  (9  mm)  notches  were  cast  into  the  center  of  12  x  12  x  1 
inches  (305  x  305  x  25  mm)  plates.  The  mortar  used  to  cast  the  specimens  had  a 
cement/sand/water  ratio  by  weight  of  1.0/2.6/0.65.  This  specimen  and  loading  type 
were  selected  because  they  provide  a  range  of  stable  crack  growth  (i.e.,  K,  decreases  as 
the  crack  length  increases).  This  enables  the  making  of  holograms  for  several  crack 
lengths  without  the  use  of  a  closed-loop  system.  This  is  important  because  the  servo- 
control  unit  would  introduce  vibrations  beyond  the  holographic  limitations. 

Three  by  six  inch  (  76  x  152  mm)  cylinders  were  also  cast  to  determine  the 
compressive  strength  and  initial  modulus  of  elasticity.  Two  cylinders  with  1/2  inch  (13 
mm)  strain  gages  were  tested  in  a  closed-loop  testing  system  under  strain  control  the  day 
after  the  center-notched  specimen  test  was  performed. 
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Figure  3  Center  Notched  Mortar  Plate  Specimen 

The  specimens  were  loaded  in  a  20-Kip  (90  KN)  screw-driven  loading  frame 
mounted  on  a  vibration  isolation  table  (Figure  4).  Two  semicylindrical  pieces  were 
attached  to  the  loading  pieces  to  apply  the  load  on  the  one  inch  diameter  loading  hole. 
The  specimens  were  loaded  in  parallel  with  four  strain-gaged  steel  bars  to  provide  extra 
stability  to  the  crack  growth. 

Holographic  interferometry  was  used  to  measure  crack  opening  displacements  and 
the  strain  field  on  a  2  inch  x  3  inch  (51  mm  x  76  mm)  area  in  front  of  each  notch  (that  is, 
top  and  bottom;  Figure  3).  It  is  important  to  notice  that  during  the  experimental  program, 
before  any  holograms  were  made,  the  specimens  were  placed  under  a  800  pound  (3.6 
KN)  preload.  This  preload  is  necessary  to  keep  the  specimen  optically  stable. 
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Figure  4  Loading  Arrangement 


C.  HOLOGRAPHIC  INTERFEROMETRY  TEST 

A  hologram  is  a  photographic  recording  that  when  illuminated  by  an  appropriate 
reference  light  reproduces  a  three  dimensional  image  of  an  object,  if  the  same  object  is 
recorded  in  two  different  states  of  deformation,  the  combined  reconstructed  image  is 
covered  with  dark  and  bright  bands  (interference  fringes)  that  represent  the  displacements 
between  the  two  states  of  deformation.  Split-beam  transmission  holograms  (Leith- 
Upatnieks  type,  Reference  12)  were  used.  A  10  mW  He-Ne  laser  with  a  wave  length  of 
24.9  /iin  (633  nm;  red)  and  coherence  length  (A)  of  approximately  3  inches  (76  mm)  was 
used  as  the  source  of  light 

Each  interference  pattern  represents  the  component  of  the  displacements  onto  the 
sensitivity  vector  k  (Reference  13).  In  the  optical  setup  shown  in  Figure  5  all  the  optical 
paths  are  on  a  plane  of  constant  height  (y  =  constant  in  Figure  3).  Therefore,  the 
sensitivity  of  the  optical  system  to  vertical  displacements  is  small.  In  addition,  the  loading 
configuration  shown  in  Figure  4  restrains  displacements  in  the  y  direction.  Thus,  two 
displacement  components  are  sufficient  to  describe  the  motion  of  any  point  on  the 
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specimen,  namely  the  in  plane  component  (x)  and  the  out  of  plane  component  (z).  To 
obtain  a  two-dimensional  displacement  field  (in  plane  and  out  of  plane),  in  principle,  two 
interferograms  are  required.  However,  since  there  is  no  sign  associated  to  the  fringe 
counts  a  third  interferogram  is  necessary.  These  requirements  are  discussed  in  more 
detail  in  a  subsequent  section. 

The  sensitivity  vector  (k),  defined  as  the  bisector  between  the  illumination  and 
observation  directions,  can  be  changed  by  either  moving  the  observation  or  the 
illumination  point.  Figure  5  is  a  schematic  view  of  the  optical  setup  with  three  illumination 
directions  (R,L  and  C)  used  for  this  study. 

For  each  load  step,  a  different  hologram  is  recorded  from  each  illumination 
direction  (holograms  R,L  and  C).  The  holograms  are  developed  and  real-time 
interferometry  is  used  to  apply  an  adequate  incremental  displacement.  One  of  the 
holograms  is  positioned  on  the  plate  holder  and  illuminated  with  the  reference  beam  while 
the  object  is  illuminated  with  the  corresponding  object  beam  (for  example,  R  object  beam 
for  hologram  R).  When  the  observer  looks  at  the  object  through  the  holographic  plate 
the  interference  between  the  real  object  and  the  holographic  image  is  viewed.  Load  is 
then  applied  to  the  specimen  until  the  fringe  pattern  reaches  the  desired  density 
(approximately  25  fringes  across  the  specimen).  A  coarse  fringe  pattern  would  result  in 
inaccurate  displacement  calculations  as  it  will  be  explained  later  and  a  fine  fringe  pattern 
would  make  the  fringe  counting  difficult  since  the  fringe  spacing  would  approach  the 
resolution  limit.  A  new  set  of  holograms  is  then  made  on  different  photographic  plates. 
When  two  holograms  corresponding  to  the  same  illumination  direction  and  consecutive 
load  steps  are  mounted  together  on  the  plate  holder  and  illuminated  with  the  reference 
beam  a  holographic  interferogram  is  reproduced.  The  above  method  is  called  sandwich 
hologram  interferometry  (Reference  14). 

The  incremental  nature  of  the  holographic  measurements,  make  it  necessary  to 
make  a  series  of  sandwich  holograms  at  continuous  load  levels.  If  a  load  increment  is 
not  recorded,  none  of  the  future  deformation  states  can  be  computed. 
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D.  IMAGE  PROCESSING 
1.  Image  Enhancement 

The  holographic  interferograms  were  acquired  into  an  image  analysis  system  using 
a  slow-scan  high-resolution  camera.  A  schematic  arrangement  for  reconstruction  of  the 
holographic  images  and  acquisition  into  the  image  analysis  system  is  shown  in  Figure  6. 
The  acquired  image  (grey  image)  contains  a  range  of  intensities  (0-255)  representing 
different  grey  levels.  The  image  is  represented  by  a  512x512  2-d  array  of  integers,  each 
corresponding  to  the  average  light  intensity  across  an  element  area  referred  to  as  pixel 
(Picture  Element).  The  size  of  a  pixel  determines  the  spatial  resolution. 

The  bright  and  dark  bands  of  the  holographic  fringes  can  be  represented  as  a 
binary  image.  In  a  binary  image  the  pixel  intensities  assume  only  two  values,  ON  (or  1) 
for  perfect  white  and  OFF  (or  0)  for  perfect  black.  Variations  of  intensity  of  the  object 
beam  throughout  the  surface  of  the  object,  make  possible  that  a  bright  band  at  a  section 
of  the  image  will  have  a  lower  intensity  magnitude  than  the  dark  bands  at  other  segments 
of  the  image.  Standard  methods  of  image  processing  such  as  thresholding  from  a  certain 
pixel  intensity  are  thus  inapplicable.  This  situation  is  represented  by  the  plot  in  Figure  7a 
which  shows  the  intensity  variations  along  a  cross  section  of  a  holographic  grey  image. 

If  I,  and  l2  are  the  intensities  of  the  individual  images  in  a  sandwich  hologram,  the 
intensity  of  the  resulting  interference  pattern  is: 

/  »  (  /,  +  4  )  +  2  (  /,  4  )1/2  cos  A<|>  00 

From  this  equation  it  is  apparent  that  the  intensity  value  of  any  pixel  can  be 
separated  into  two  terms.  One  is  the  average  intensity  of  the  image  (I,  +  y,  and  the 
second  depends  on  the  phase  difference  (A<p)  between  the  two  wave  fronts.  The 
separation  of  this  two  terms  can  be  achieved  through  the  use  of  the  image  analysis 
system. 

Assume  that  the  grey  image  of  the  sandwich  holography  interferogram  is 
represented  as  IMAGE_1  where  as  IMAGE_2  is  the  average  of  the  independent  wave 
fronts  acquired  separately  (Figure.  7a).  Subtraction  of  IMAGE_2  from  IMAGE_1  results 


in  the  isolation  of  the  interferometric  effect,  represented  as  IMAGE_3  in  Figure  7b.  Since 
negative  values  are  not  allowed  as  meaningful  intensity  values,  IMAGE_3  results  in  zero 
for  all  the  pixels  on  dark  bands.  The  bright  bands  retain  a  positive  magnitude.  By 
selecting  all  positive  values  a  binary  image  is  created  (Figure  7c). 


2.  Fong®  Count 

1MAGE_3  is  the  enhanced  interferometric  image  in  which  the  intensity  value  is  zero 
for  dark  bands  and  one  for  bright  bands.  Further  image  processing  allows  to  assign 
intensity  values  to  each  fringe  corresponding  to  the  fringe  order. 

A  pixel  P  is  selected  on  the  binary  image  (IMAGE_3).  A  flooding  procedure  which 
assigns  an  arbitrarily  selected  intensity  value  (Y)  to  all  the  pixels  connected  to  Pixel  P  is 
performed.  Thus,  ail  the  pixels  on  the  band  containing  P  are  assigned  the  intensity  value 
Y.  The  same  procedure  is  repeated  for  consecutive  bands.  The  intensity  value  Y  is 
incremented  by  one  for  each  band. 

Even  though  this  process  has  been  greatly  automated,  user  interaction  is 
necessary  for  the  definition  of  the  fringe  count  direction.  In  Figure  7c  line  A-B  indicates 
the  direction  of  positive  fringe  count  (e.g.,  the  fringe  order,  or  intensity  value,  increases 
from  A  to  B).  The  above  flooding  process  is  repeated  until  the  whole  area  of  interest  has 
been  included.  In  Figure  8  flooding  trajectories  A-B,  C-D  and  E-F  have  to  be  followed. 

Figure  8  shows  a  fully  processed  interferogram  in  which  same  intensity  values 
correspond  to  same  fringe  orders.  Note  that  there  is  a  different  intensity  level  for  each 
half  a  fringe  (that  is,  different  intensity  for  dark  and  bright  bands).  The  fringe  count 
between  any  two  points  can  now  be  easily  obtained  by  computing  half  the  difference  of 
the  corresponding  intensity  levels  (e.g.,  fringe  count  between  A  and  B  in  Figure  8  is  (90  - 
65)/2  =  12.5).  For  further  refinement  on  the  fringe  count,  the  fringe  order  of  the  point 
A  in  Figure  8  is  computed  according  to  Equation  (2)  (see  inset  of  Figure  8). 

NP*  (2) 

2  2  A/P 
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Figure  8  Fringe  Count  by  Difference  in  Intensity  Value 
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Even  though  the  use  of  the  image  analysis  system  results  in  very  fast  fringe  counts, 
speed  is  not  the  main  advantage  of  this  technique.  The  isolation  of  the  interferometric 
effect  on  the  holographic  images  and  the  binarization  process  result  in  accurate  and 
objective  fringe  counts.  The  gradual  change  in  intensity  from  a  dark  to  a  bright  band 
make  it  difficult  to  obtain  consistent  fringe  counts  when  conventional  (by  hand)  methods 
are  used. 


E.  EVALUATION!  OF  HOLOGRAMS 

The  holodiagram  (Reference  1 5)  shown  in  Figure  9a  is  a  simple  geometrical  tool 
for  computation  of  displacements  based  on  interferometric  readings.  In  Figure  9a  the 
distance  i-P-0  from  illumination  point  to  object  to  observation  point  (usually  taken  as  the 
position  of  the  holographic  plate)  remains  constant  if  the  object  (P)  moves  along  ellipse 
E1  (or  an  ellipsoid  in  the  most  general  three  dimensional  case)  with  foci  at  illumination  and 
observation  points.  Now,  define  ellipse  E2  such  that  the  distance  l-P-0  increases  by  X/2 
(180°  phase  difference).  An  interferogram  of  the  object  that  moved  from  any  point  on 
ellipse  E1  to  any  point  on  ellipse  E2  will  show  a  dark  band  (half  a  fringe)  due  to  destructive 
interference  of  the  wave  fronts.  In  general,  the  number  of  fringes  is  proportional  to  the 
increase  (or  decrease)  of  the  distance  I-P-0  (Equation  (3)). 

A  (  l-P-0  )=n  X  (3) 

where:  A  (l-P-O)  increment  of  distance  l-P-0 

n  fringe  count  (  one  fringe  =  one  dark  band  and  one  bright  band) 


E  Ellipce  for  constant  l-P-0 

Ej+  Ellipce  for  +1/2  fringe  (  A-P-0  =  +  \/2) 

Es  Ellipce  for  -1  /2  fringe  (  zl-P-0  =  -  V2) 

I  Illumination  Point 
O  Observation  Point 
P  Point  on  the  Object 

Several  possible  displacement 
vectors  for  1/2  fringe. 


Figure  9  Evaluation  of  Interferograms 

a)  Holodiagram 

b)  Displacement  Vector  as  the  Intersection  of  Two  Holodiagrams 
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is  the  cumulative  of  the  CODs  for  all  the  previous  load  levels.  Two  typical  crack  profiles 
included  in  Figure  10  show  that  there  is  a  good  match  in  the  crack  profiles  obtained  with 
both  holograms. 

An  alternative  solution  is  to  assume  a  two  dimensional  displacement  solution.  The 
displacement  is  computed  as  the  difference  in  position  of  the  intersection  of  the  ellipses 
from  two  illumination  directions  (Figure  9b).  The  holographic  fringe  counts  can  be 
combined  in  three  possible  ways  (  R  and  L,  R  and  C,  C  and  L).  As  mentioned  before, 
the  third  fringe  count  is  necessary  to  select  the  correct  value  among  the  four  possible 
displacement  solutions  for  each  combination.  Figure  10  shows  a  crack  profile  computed 
using  all  three  combinations.  There  is  an  improvement  in  the  scatter  on  the  COD 
calculations  with  respect  to  the  Nelson  and  McCrickerd  results.  Also,  note  that 
considering  the  out  of  plane  contribution  to  the  fringe  counts  results  in  smaller  CODs, 
even  for  an  experimental  setup  designed  for  Mode  I  crack  propagation.  This  means  that 
there  exists  a  measurable  component  of  the  crack  face  displacement  exists  in  the  out-of¬ 
plane  direction  (Mode  III)  which  can  explain  the  presence  of  compressive  strain  behind 
the  crack  tip  as  reported  by  Bruhwiler.  The  dotted  line  in  Figure  10  represents  the  crack 
sliding  displacement  (CSD)  in  the  XZ  plane. 

2.  Strain  Fields. 

Fringe  counts  relative  to  the  left  face  of  the  crack,  or  the  symmetry  line  for  point 
beyond  the  crack  tip,  were  taken  at  every  point  on  a  3  inch  x  3  inch  (76  mm  x  76  mm) 
grid  (1 .5  inches  on  each  side  of  the  crack  line  and  3  inches  in  front  of  the  notch  tip)  every 
1  /8  inches  (3  mm)  for  all  three  holographic  images  (R,L  and  C).  A  computer  program 
was  used  to  calculate  the  relative  displacement  of  every  point  relative  to  the  left  face  of 
the  crack  by  locating  the  intersection  point  of  the  holodiagrams  corresponding  to  each 
illumination  direction.  A  third-order  best-fit  polynomial  was  obtained  for  the 
displacements  along  lines  perpendicular  to  the  crack  line.  For  the  sections  crossing  the 
crack  a  different  polynomial  was  obtained  for  the  left  and  right  side  of  the  crack.  The 
strain  field  was  obtained  by  differentiation  and  was  evaluated  from  x=-1  in  to  x=  +  1  in 
every  1  /8  in  (3  mm)  at  every  cross  section,  where  x  is  the  distance  from  the  crack  line 
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Figure  10  Holographically  Measured  Crack  Profiles 


and  is  positive  to  the  right.  For  the  sake  of  brevity,  results  for  one  of  the  cracks  on  one 
of  the  specimens  are  shown.  These  results  can  be  considered  representative  of  the  six 
cracks  examined.  Figure  11  shows  the  strain  field  (e^)  around  the  bottom  crack  for 
crack  lengths  of  1 .375  in  (35  mm),  1.875  in  (48  mm)  and  2.25  in  (57  mm),  corresponding 
to  1324  lb  (6.0  KN),  1461  lb  (6.6  KN)  and  1601  lb  (7.2  KN)  of  applied  load.  The  strain 
fields  shown  in  Figure  11  include  only  the  measured  additional  strain  after  the  initial 
preload  state. 

In  general,  a  high  strain  (>100jue)  region  can  be  found  around  and  in  front  of  the 
crack  tip  with  a  maximum  value  of  approximately  300  ne.  Behind  the  crack  tip  there  is  a 
zone  of  negative  strain  which  suggests  relaxation  of  tensile  strain  with  respect  to  the  initial 
preload  stage. 

Figure  12  shows  the  strain  distribution  along  sections  parallel  to  the  crack  line  for 
various  load  levels.  In  this  Figure,  a  linear  elastic  component  corresponding  to  the  initial 
preload  has  been  added  in  an  attempt  to  provide  a  better  picture  of  the  total  state  of 
strain.  Similar  observations  as  those  made  by  Cedolin,  Dei  Poli  and  lori  (Reference  7) 
can  be  made.  Figure  1 2  shows  that  the  material  on  both  sides  of  the  crack  unloads.  The 
strain  concentrations  near  the  crack  tip  and  the  unloading  are  greater  on  sections  closer 
to  the  crack  line.  Although  there  is  unloading  behind  the  crack  tip,  the  strain  level  does 
not  immediately  reach  zero;  this  suggests  the  existence  of  tension  being  transmitted 
through  the  crack. 

F.  DEFINITION  OF  THE  FRACTURE  PROCESS  ZONE 

Cedolin,  et  al.  (Reference  7)  suggested  the  use  of  isostrain  contours  to  enclose 
zones  of  microcracking  and  he  referred  that  zone  as  Fracture  Process  Zone  (FPZ).  He 
assumed  that  the  region  where  the  strain  value  exceeded  150  ne  was  FPZ.  In  this  study 
it  was  assumed  that  the  strain  value  of  100  Ate  is  a  more  reasonable  limit  to  define 
nonlinearity.  The  zone  where  the  measured  strain  values  exceed  100  tie  are  shown  in 
Figure  13  for  three  crack  lengths.  Also  shown  in  the  Figure,  the  shaded  areas  are  the 
zones  behind  the  crack  tip  where  tensile  strain  was  detected.  The  end  of  these  zones 
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Distribution  Along  Section  Parallel  to  the  Crack  Line 


correspond  to  the  tensile  strain  value  of  zero  (see  Figure  11).  There  appears  to  be  a 
traction  free  zone  of  approximately  0.25  inch  in  front  of  the  original  notch  tip  regardless 
of  the  crack  length  and  COD  indicated  by  wc  in  Figure  13).  This  agrees  with  the  model 
suggested  by  Cook,  et  al.  (Reference  7)  where  it  is  postulated  that  cohesive  forces  are 


strain  drops  from  500  jue  at  the  crack  tip  to  zero  within  1  /8  in  behind  the  crack  tip.  In 
contrast,  the  observed  strain  relaxation  shown  on  Figures  11  and  12  follows  a  gradual 
descent. 

In  Figure  15,  the  differences  in  strain  fields  between  the  LEFM  solution  and  the 
experimentally  measured  values  are  shown.  It  is  arbitrarily  assumed  that  the  differences 

t 

the  tensile  strain  predicted  by  LEFM  was  higher  while  the  negative  values  mean  that  th 


h 


marked  A  in  Figure  15  are  the  regions  where  the  LEFM  solution  predicts  tensile  strain 
values  of  60  jue  or  more  higher  than  those  measured  holographically.  Note  that  Zone  A 
remains  essentially  constant  regardless  of  crack  length.  Since  this  zone  is  relatively  small 
and  does  not  change  with  crack  length  a  model  based  on  modified  LEFM  may  be  a 
possible  approach.  This  is  discussed  in  the  next  section.  Zone  B,  with  negative  values, 
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Figure  13  Fracture  Process  Zone 
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typically  found  behind  the  crack  tip,  shows  that  the  relaxation  of  the  observed  strain  is  60 
Me  less  than  that  in  the  elastic  solution.  Zone  B  can  be  defined  as  the  wake  of  the 
fracture  process  zone  (WFPZ)  and  encloses  the  area  where  extensive  microcracking  has 
been  developed  and  tensile  forces  are  still  transmitted  through  the  crack.  Behind  this 
zone,  the  fracture  process  approaches  a  traction-free  condition.  Similar  conclusions 
regarding  the  growth  and  propagation  of  the  FPZ  were  reached  analytically  by  Bazant  and 
Gettu  (Reference  17)  if  the  crack  tip  is  assumed  behind  the  WFPZ. 


G.  CLOSING  PRESSURE  ANALYSIS 

A  cohesive  type  of  model  has  been  used  to  simulate  effect  of  fracture  process 
zone.  A  closing  pressure  is  applied  to  the  faces  of  the  crack  and  it  is  assumed  that  the 
closing  pressure  is  a  function  only  of  crack  opening  displacement.  To  determine  if  such 
a  model  can  accurately  predict  the  measured  strain  field,  the  same  finite-element  program 
with  singular  elements  was  used.  A  bilinear  closing  pressure  vs.  COD  relationship  was 
used.  This  relationship  was  shown  by  Miller  (Reference  18)  to  reasonably  accurately 
predict  the  holographically  measured  crack  profiles  for  similar  specimens  made  with 
identical  mortar  mix  proportions.  In  general,  an  iteration  procedure  is  necessary  to 
perform  the  analysis  based  on  a  closing  pressure  vs.  COD  relation.  However,  since  the 
CODs  are  expected  to  converge  close  to  the  holographically  obtained  values,  these  can 
be  used  to  compute  the  closing  pressure.  Thus,  the  iteration  process  is  not  necessary. 

Figure  16  shows  a  comparison  of  the  holographically  observed  crack  profiles,  the  LEFM 
solution  and  the  predicted  crack  profiles  using  the  bilinear  closing  pressure.  The 
predicted  crack  profiles  in  Figure  16  show  a  good  correlation  with  the  experimental 
values. 

Figure  1 7  shows  the  difference  between  the  strain  field  predicted  using  the  closing 
pressure  approach  and  the  experimentally  observed  values.  There  is  a  good  correlation 
of  the  results  and  the  zones  with  significant  deviation  from  the  experimental  results  are 
practically  eliminated  (±  60  j ue).  Figure  18  reports  the  computed  values  of  the  resultant 
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CRACK  LENGTH. 
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through  the  crack  growth.  Thus,  a  constant  K,  =  Klc  type  of  criteria  can  be  used  to 
predict  crack  growth. 


h.  eoweiusiGMis 

1.  Sandwich  hologram  interferometry  with  multiple  sensitivity  vectors  can  be  used  to 
measure  crack  opening  displacements  and  strain  fields  in  mortar  specimens  under 


opening  displacements  the  effect  of  out  of  plane  motion  should  be  taken  into  account. 

2.  Digital  image  analysis  facility  provided  a  faster,  more  consistent  and  more  accurate 
method  of  fringe  count  than  the  manual  method. 

3.  There  is  a  region  of  high  tensile  strain  that  moves  with  the  tip  of  the  propagating 
crack.  The  material  behind  the  crack  tip  experiences  strain  relaxation.  However,  the 


5.  A  bilinear  closing  pressure  vs.  COD  cohesive  crack  model  has  been  proposed 
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Difference  Between  Holographic  and  Bilinear  Closing  Pressure  Si 
(a)  Load  =  1324  lb.  (b)  Load  =  1461  lb.  (c)  Load  =  1601  lb. 


CRACK  LENGTH,  a  (mm) 


Figure  18  Apparent  Mode  I  Stress  Intensity  Factor(Kj)  vs.  Crack  Length  (a) 


SECTION  HI 


CRACK  INITIATION  AND  PROPAGATION  IN  DISK  SPECIMEN 

(BRAZILIAN  TEST) 


A.  INTRODUCTION 

A  common  method  for  estimating  the  tensile  strength  of  concrete  is  the  splitting  test. 
The  test  is  carried  out  by  loading  a  cylinder  on  its  side  in  diametral  compression.  This 
indirect  method  to  determine  tensile  strength  is  attributed  to  Carneiro  (Reference  19) 
and  was  adopted  the  national  standard  specifications  of  Brazil,  thus  the  name  "Brazilian 
Test."  According  to  the  theory  of  elasticity,  this  loading  configuration  produces  a  nearly 
constant  stress  along  approximately  75  percent  of  the  diameter  (References  20  and 
21).  Although  the  state  of  stress  is  not  uniaxial  (significant  compressive  stress  in  the 
transverse  direction)  and  there  is  redistribution  of  stress  due  to  nonlinearity  of  concrete, 
the  tensile  stress  value  at  failure  has  proven  to  be  a  useful  measure  of  the  tensile 
strength.  Figure  19  shows  the  failure  sequence  in  the  disk  specimen  (Reference  22). 
The  purpose  of  this  study  is  to  evaluate  the  distribution  of  deformations  on  the  surface 
of  a  disk  specimen  as  it  is  loaded  and  determine  the  failure  sequence. 

B.  TESTING  PROCEDURE 

1.  Mechanical  Configuration 
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frame  shown  in  Figure  20.  The  use  of  such  a  large  loading  frame  was  selected  because 
it  includes  a  heavy  platform  which  was  used  as  a  vibration  free  table  for  the  holographic 


attempts  were  made  using  cardboard  strips  to  distribute  the  load  over  the  desired  area, 
but  the  deformation  of  the  cardboard  during  the  application  of  the  load  was  orders  of 
magnitude  higher  than  the  deformation  of  the  specimen,  resulting  in  holographic 
recording  reflection  the  deformation  of  the  loading  strip  instead  of  the  specimen.  The 
specimens  were  placed  between  the  loading  platens  using  steel  spacers  designed  to 
support  the  instrumentation  and  keep  the  specimen  12  inches  (309  mm)  above  the 
surface  of  the  loading  frame  platform  (Figure  22). 


Figure  20  120  Kip  Capacity  Loading  Frame 


The  compressive  strength  and  modulus  of  elasticity  of  each  batch  of  mortar  were 
determined  by  testing  in  uniaxial  compression  2-inch  x  4-inch  (50.8mm  x  101.6mm) 
cylinders  at  comparable  age.  The  same  120-Kip  loading  frame  was  used.  The  dual  LVDT 
transducer  shown  in  Figure  23  was  used  to  measure  axial  displacement  and  the  applied 
load  was  measured  using  the  AP  transducer  of  the  loading  frame.  The  specimens  were 
tested  at  constant  rate  of  deformation  such  that  peak  load  was  reached  in  about  20 
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during  testing.  The  AP  transducer  existing  in  the  loading  frame  was  not  adequate  to 
measure  load  levels  in  the  order  of  1  to  2  Kips.  The  cylindrical  load  cell  shown  in  Figure 
22  was  designed  and  fabricated  for  accurate  load  measurement  and  holographic  stability. 
The  load  cell  was  designed  for  a  maximum  load  of  6000  lb  in  compression.  A  full 
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Figure  24  Stress-Strain  Curves  for  Uniaxial  Compression  Cylinders 


Wheatstone  bridge  was  used  on  the  load  cell  to  measure  axial  load  only  (eliminate 
eccentricity  effect)  and  compensate  for  temperature. 

The  load  point  displacement  was  measured  using  the  same  dual  LVDT  transducer 
(designed  and  fabricated  during  the  course  of  this  investigation)  used  to  obtain  the 
modulus  of  elasticity.  As  shown  in  Figure  23,  the  dual  LVDT  transducer  was  attached  to 
the  spacer  blocks. 
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AP  and  stroke  readings  were  recorded. 


3.  Control  Signal 

Preliminary  tests  were  performed  in  order  to  determine  the  most  suitable  control 
system.  Figure  26  shows  a  schematic  of  the  results.  In  the  load-load  point  displacement 
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Figure  25  Diametral  Displacement  Transducer 


plot  there  are  three  zones:  (Zone  A)  linear  increase  of  the  load,  (Zone  B)  transition  stage 
where  the  first  crack  propagates  and  (Zone  C)  plateau  where  secondary  cracking  and 
final  failure  occurs.  The  main  drawback  of  the  displacement  control  system  is  that  the 
transition  stage  is  unstable  (drop  of  the  load  at  constant  displacement). 

The  load-DD  curve  also  shown  in  Figure  26  includes  the  same  three  zones.  The  DD 
in  zone  A,  corresponding  to  the  elastic  stage,  is  almost  zero.  If  this  signal  was  used  for 
control,  a  sudden  increase  in  the  load  would  occur  prior  to  cracking.  Zones  B  and  C 
show  the  desired  stable  behavior  (i.e.  continuous  increase  of  the  DD)  after  cracking  has 
started. 

A  new  control  system  was  designed  in  order  to  have  stable  behavior  throughout  the 
loading  sequence.  Figure  26  shows  that  in  Zone  A  the  load  increases  while  the  DD 
remains  constant  while  in  zones  B  and  C  the  load  variation  is  small  and  the  DD  increases. 
By  selecting  a  linear  combination  of  the  load  and  DD  signals  the  desired  response  is 
obtained  as  illustrated  in  Figure  26.  Note  that  the  load  +  DD  signal  continuously  increase 
in  all  stages  of  the  loading  sequence. 

Since  the  addition  of  two  signals  cannot  directly  be  used  for  control  by  our  analog 
closed  loop  control  system,  a  feedback  control  system  was  designed.  Figure  27  shows 
a  block  diagram  of  the  testing  system.  In  this  diagram  the  shaded  areas  correspond  to 
the  components  developed  during  the  course  of  this  investigation.  Figure  28  shows  a 
block  diagram  of  the  load  +  DD  feedback  control  amplifier.  The  main  components  of  the 
control  amplifier  are  two  channels  of  preamplifiers,  a  signal  summator  and  three  output 
drivers.  The  preamplifiers  are  used  to  select  the  relative  weight  of  the  load  and  DD 
signals.  The  summator  is  used  to  add  the  preamplified  signals  and  the  output  drivers  are 
used  to  control  the  calibration  factors  and  provide  separate  load,  DD  and  load  +  DD 
signals  for  monitoring  and  recording. 


C.  HOLOGRAPHIC  IWTERFEROf 


TEST 
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Figure  26  Control  Signal 
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Figure  28  Load+DD  Feedback  Control  Amplifier 


To  measure  three  dimensional  displacement  fields  a  three  sensitivity  vector  holographic 
arrangement  was  required  (see  Section  III.E).  Split  beam  transmission  sandwich 
holography  (see  Section  II. C)  was  used.  The  complete  optical  arrangement  is  shown  in 
Figures  29  and  30.  For  better  use  of  the  limited  space  two  horizontal  planes  were  used 


Figure  30  Optical  Arrangement  for  Three  Sensitivity  Vectors 


for  the  optical  arrangement  (as  opposed  to  the  more  conventional  single-plane  setup). 
For  easier  understanding  of  the  operation  of  the  holographic  setup,  Figure  31  shows  only 
the  most  significant  components  of  the  arrangement.  To  obtain  three  different  sensitivity 
vectors,  a  common  reference  (r)  beam  and  three  different  object  beams  (R,L,D)  were 
used. 

As  noted  in  Section  II. C,  holographic  interferometry  experiments  must  be  performed 
under  very  stable  conditions.  Conventionally,  noise,  light,  air  currents  and  vibrations  are 
isolated  from  the  system  through  the  use  of  devices  such  as  vibration  isolation  tables  and 
noise  dampers.  However,  during  this  investigation,  the  use  of  the  servo-controlled  dosed 
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Figure  29  Optical  Arrangement  for  Three  Sensitivity  Vectors  (Schematic) 


47 


loop  testing  system  was  essential  for  the  appropriate  mechanical  behavior.  Extreme  care 
was  placed  in  the  adaptation  of  the  laboratory  for  holographic  recordings.  The  hydraulic 
pump  was  placed  on  rubber  dampers  to  reduce  vibrations  transmitted  through  the  floor 
to  the  loading  frame.  A  tent  made  of  light  tight  fabric  was  built  around  the  testing  frame 
to  avoid  air  currents  and  provide  dark  room  conditions.  The  plate  holder  shown  in  Figure 
32  was  designed  to  record  three  holograms  on  the  same  photosensitive  plate  which 
eliminates  the  vibrations  induced  by  changing  plates  after  each  exposure.  Also,  after  the 
plate  has  been  developed  and  the  real-time  technique  is  used  to  determine  the  size  of  the 
load  increment,  all  three  holograms  can  be  viewed  at  the  same  time.  Thus,  there  is  no 
need  to  “guess"  the  critical  component  of  the  displacements.  A  Richo  XR-M  camera 
without  lens  was  programmed  to  make  three  exposures  at  30-second  intervals  and  placed 
in  front  of  the  laser  head  to  act  as  a  shutter.  For  each  exposure,  only  one  of  the  object 
beams  was  active  (the  other  two  were  blocked  with  cardboard  screens).  In  addition,  all 
holographic  experiments  were  performed  overnight  or  on  weekends  to  avoid  the  noise 
generated  by  normal  use  of  the  building.  Before  an  experiment  was  performed,  all  the 
mechanical  components  were  tightened,  the  transducers  were  mounted  and  the  electronic 
systems  were  operated  for  24  hours  to  allow  for  all  the  components  to  "seat“  into  a  stable 
position.  The  electronic  systems  were  started  24  hours  in  advance  and  the  hydraulic 
system  was  started  3  hours  in  advance  in  to  reach  thermal  stability  of  the  system.  Since 
the  duration  of  a  holographic  experiment  is  approximately  8  hours  and  holographic 
stability  must  maintained  throughout  the  whole  process,  even  with  all  these  precautions, 
the  success  ratio  of  the  holographic  experiments  was  approximately  1/10. 

Two  reference  objects  were  included  in  the  holographic  recordings  (Figure  22).  One 
was  glued  to  the  spacer  block  and  the  other  to  the  specimen.  The  reference  object  on 
the  spacer  block  was  used  to  determine,  in  case  of  low  quality  holograms,  if  the  source 
of  motion  was  due  to  variable  load  applied  on  the  specimen  or  vibration  on  the  loading 
platform  itself.  The  second  reference  object,  glued  to  the  specimen,  was  used  to  evaluate 
the  rigid  body  displacement. 
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D.  IMAGE  PROCESSING  (ELIMINATION  OF  RIGID  BODY  MOTION) 

Holographic  interferometry  is  sensitive  to  displacements,  including  both  rigid  body 
motion  and  displacements  induced  through  deformation.  In  the  study  presented  in  Section 
II,  the  effect  of  the  rigid  body  motion  was  neglected  because  the  mechanical  configuration 
allowed  only  rotation  about  the  axis  of  loading.  If  the  horizontal  component  of  the  relative 
displacement  between  Points  A  and  B  in  Figure  33  is  calculated,  the  rigid  body  motion 
displacement  of  Point  A  is  the  same  as  that  of  Point  B.  Thus,  an  arbitrary  point  of  zero 
displacement  can  be  selected  and  incremental  displacements  with  respect  to  this  point 
suffice  to  compute  the  deformation  field. 

However,  in  this  study,  three-dimensional  deformation  fields  are  computed.  Unless 
the  point  of  zero  displacement  (with  respect  to  the  photographic  plate)  is  on  a  known 
position  on  the  surface  of  the  object,  elimination  of  rigid  body  motion  is  necessary. 

Abramson  (Reference  13)  showed  that  the  rigid  body  motion  can  be  eliminated  from 
an  interferogram  by  subtracting  the  fringe  pattern  corresponding  to  the  rigid  body  motion 
alone.  The  rigid  body  motion  can  also  be  eliminated  by  rotating  the  sandwich  holograms. 
However,  as  also  noted  by  Abramson,  the  rigid  body  motion  component  provides 
additional  information  regarding  the  direction  of  the  displacements. 

The  surface  of  the  mortar  specimen  was  subjected  to  both  deformation  and  rigid  body 
motion  between  holographic  exposures,  but  the  reference  object  glued  to  the  side  of  the 
specimen  was  only  subjected  to  rigid  body  motion.  Figure  34  shows  an  original  gray 
image  of  an  interferogram  of  the  disk  specimen.  Note  that  both  the  reference  objects 
contain  holographic  fringes  (rigid  body  motion).  Also  note  that  the  fringe  patterns  on  both 
reference  objects  are  different  which  suggests  that  there  was  displacement  of  the 
specimen  with  respect  to  the  loading  block.  The  fringe  orders  were  evaluated  as  in 
described  in  Section  II. D. 2  resulting  in  Figure  35.  In  this  image,  the  fringe  orders  have 
integer  values  (fringe  numbers).  Note  that  the  reference  object  was  also  colorized.  The 
fringe  order  between  the  reference  object  and  the  mortar  specimen  does  not  have  to  be 
continuous  but  the  direction  of  positive  fringe  count  has  to  be  consistent. 
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Figure  34  Holographic  Interferogram  of  Disk  Specimen  (Original 
Image).  Hologram  3 


Direct  application  of  the  Moire  pattern  method  (superimpose  a  grid  of  the  same 
orientation  as  density  as  that  on  the  reference  object  on  the  entire  image)  results  on 
elimination  of  the  rigid  body  motion  but  the  transition  between  bright  and  dark  fringes 
becomes  ambiguous  (Figure  36).  A  computer  program  was  written  to  mathematically 
subtract  the  fringe  patterns  based  on  the  orientation  and  density  of  the  reference  grid  and 
the  distance  to  the  point  of  zero  displacement  (fixture  between  specimen  and  reference 
object).  The  fringe  order  at  any  point  P  on  the  object  after  fringe  superposition  is  given 
by  Equation  (5). 

/=/,-[(  of  x  A  )•  k)  df  (5) 

where:  /  Fringe  order  after  rigid  body  motion  elimination 

/,  Original  fringe  order 

of  Unit  vector  defining  the  orientation  of  the  reference  grid 
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Figure  35  Colorized  Fringe  Pattern  Before  Elimination  of  Rigid 
Body  Motion.  Hologram  3 


df  Fringe  density  (fringes/unit  distance) 

A  Vector  from  point  of  zero  displacement  to  point  P 
k  Unit  vector  in  the  z  direction  (out  of  plane) 


Note  that  in  this  equation  the  original  fringe  orders  are  real  values  obtained  by 
Equation  (2).  Therefore,  the  fringe  orders  in  the  resulting  image  shown  in  Figure  37  are 
also  real  numbers.  This  fringe  orders  can  directly  been  used  for  displacement  calculation. 


E„  EVALUATION  OF  HOLOGRAMS  (STRAIN  FIELDS) 

In  a  holographic  setup,  if  the  distances  between  observation  points  and  illumination 
point  to  the  object  are  large  compared  to  the  magnitude  of  the  displacements,  the 
following  linear  approximation  can  be  made  to  the  exact  holodiagram  solution  (Figure  9). 
The  displacement  vector  d  can  be  decomposed  into  three  orthogonal  components,  one 
along  the  sensitivity  vector  k  (defined  as  the  bisector  of  the  angle  IPO  in  Figure  38)  and 
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Figure  36  Elimination  of  Rigid  Body  Motion  by  Direct  Moire  Patter 
Superposition.  Hologram  3 


two  perpendicular  to  k.  For  small  displacements,  motion  of  the  point  P  perpendicular  to 
k  cause  no  change  in  in  the  path  length  l-P-0  (the  perpendicular  motions  are  along  the 
ellipse  having  I  and  O  as  the  foci).  The  linear  approximation  of  the  displacement  d  is 
given  by: 

d  -  k  =  1  =  H  ~  (6) 

2  cos  a  2  cosa 

To  separate  the  tree  components  of  the  displacement  in  the  most  general  case, 
independent  holographic  readings  have  to  be  made  with  three  different  (non-coplanar) 
sensitivity  vectors  and  the  above  equation  takes  the  form: 
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Figure  37  Colorized  Image  After  Mathematical  Elimination  of  Rigid 
Body  Motion.  Hologram  3. 


In  this  equation,  the  displacement  vector  d  is  a  function  only  of  the  fringe  counts  nif 
and  the  known  wave  length  A.  and  geometrical  constants  kj  and  otj. 

A  computer  program  was  written  to  compute  the  displacements  using  the  linear 


holographic  system.  Displacements  were  evaluated  every  tenth  of  a  inch  on  the  surface 
of  the  specimen.  The  procedure  described  in  Section  li.E.2.  was  used  for  strain 
computations. 

The  load  history  is  shown  in  Figure  39.  The  plateaus  on  this  curve  correspond  to  the 
time  necessary  to  process  the  holographic  plates  during  which  the  load  was  kept 
constant.  In  this  figure,  the  stars  correspond  to  the  state  when  a  hologram  was  made 
as  shown  in  Table  2.  Figures  40  to  44  show  the  strain  fields  (e^)  at  three  stages  of 
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Figure  38  Linear  Approximation  for  Evaluation  of  Holograms 


loading.  The  linear  elastic  solution  was  obtained  using  FRANC  (FRacture  ANalysis  Code), 
an  interactive  finite-element  program  that  models  crack  tips  though  six-node  quarter  point 
node  traingular  elements  and  provides  automatic  remeshing  and  crack  propagation 
(Reference  23).  Note  that  both  the  holographic  and  the  linear  elastic  strain  fields 
correspond  to  the  increment  of  strain  from  the  initial  preload  state  of  383  pounds  to  the 
given  load  level.  Holograms  3  (Figures  37  and  40)  give  a  typical  strain  distribution  in  the 
"elastic"  range.  The  measured  values  show  the  same  trends  and  maximum  tensile  strain 
of  100  /i€  as  the  numerical  solution.  However,  the  holographic  results  show  a 
concentration  of  strain  near  the  point  of  application  of  the  load  at  the  bottom  of  the 
specimen,  which  may  suggest  that  even  at  this  load  level  (first  crack  at  2,738  pounds) 
the  crack  initiation  location  can  be  predicted.  Holograms  11  (Figures  41  and  42) 
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TABLE  2  HOLOGRAPHIC  RECORDINGS.  SPECIMEN  HOLO-UN-1 


Holographic 

Recording 

Time 

Load 

s 

lb 

N 

P  (Preload) 

427 

383 

1,724 

1 

2,500 

648 

2,916 

2 

4,163 

913 

4,109 

3 

6,250 

1,089 

4,901 

4 

7,800 

1,281 

5,765 

5 

9,310 

1,472 

6,624 

6 

11,116 

1,752 

7,884 

7 

12,785 

1,958 

8,811 

8 

14,584 

2,179 

9,806 

9 

18,000 

2,282 

10,269 

10 

19,560 

2,517 

11,327 

11 

21,113 

2,738 

12,321 

12 

22,600 

2,959 

13,316 

13 

24,400 

3,180 

14,310 

14 

25,200 

3,533 

15,899 

15 

25,860 

3,901 

17,555 

solution.  A  second  strain  concentration  is  observed  near  the  top  of  the  specimen 
indicating  probable  crack  development  after  further  loading.  A  significant  difference 
between  this  study  and  the  center-notched  plate  specimen  presented  in  Section  II  is  that 
the  strain  levels  are  significantly  higher  for  the  cylinder  specimen  as  compared  with  the 
center  loaded  plate  specimen  even  on  the  uncracked  region.  This  can  be  attributed  to 
the  presence  of  high  compressive  stresses  in  the  transverse  direction  but  it  has  not  yet 
been  investigated.  Although  further  analysis  is  necessary,  a  preliminary  investigation 
based  on  the  elastic  solution  shows  that  the  mode  II  (shear)  stress  intensity  factor  (K,,) 
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Figure  41  Interferometric  Fringes.  Hologram  11. 


is  approximately  10  times  higher  than  K,. 

Holograms  14  (Figures  43  and  44)  were  taken  at  the  last  stable  load  level  and  can  be 
considered  the  failure  stage.  At  this  load  stage,  strain  could  not  be  obtained 
experimentally  near  the  points  of  application  of  the  load.  The  size  of  the  interferometric 
fringes  was  beyond  the  resolution  of  the  optical  system  used.  It  was  inconclusive  if  the 
crack  propagated  to  the  edges  of  the  specimen  or  not.  At  this  stage  the  crack  has  grown 
0.2  inch  towards  the  center  of  the  specimen.  For  this  crack  length,  the  elastic  solution 
shows  the  sharp  strain  relaxation  behind  the  crack  tip  typical  of  noncohesive  cracks  while 
the  experimental  results  still  show  significant  strain  behind  the  crack  tip.  At  the  crack  tip 
near  the  center  of  the  specimen,  the  LEFM  K„/K(  ratio  is  approximately  1. 

Additionally,  Figure  45  shows  two  interferograms  corresponding  to  exposures  at  3901 
pounds  of  applied  load.  Note  that  the  crack  length  is  different  for  interferograms  14D15 
and  14R15  (Down  and  Right  object  beams  respectively),  which  indicates  crack 
propagation  between  exposures  for  the  same  applied  load.  The  third  hologram  (14L15) 
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Figure  43  Interferometric  Fringes.  Hologram  14. 

was  not  recorded  because  the  crack  propagated  and  the  specimen  failed  during  its 
exposure.  The  specimen  was  loaded  beyond  the  range  of  the  displacement  contorl 
transducers  under  constant  rate  of  load  point  displacement  (stroke)  in  order  to  determine 
the  final  failure  mode.  The  observed  failure  sequence  is  shown  in  Figure  46.  A  very 
similar  failure  sequence  at  high  strain  rates  have  been  reported  by  Tedesco  et  al. 
(Reference  24) 

F.  CONCLUSIONS 

1.  Laser  holographic  interferometry  can  be  used  to  measure  three  dimensional 
displacement  fields  under  closed  loop  testing  conditions. 

2.  A  mechanical  arrangement  was  developed  to  produce  stable  crack  growth  in  disk 
specimens  (Brazilian  Test). 

3.  Image  processing  can  be  effectively  used  for  elimination  of  rigid  body  motion  from 
holographic  recordings  resulting  in  well  defined  and  accurate  fringe  orders. 
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Figure  45  Holographic  Interferograms  at  Failure  State. 
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4.  Examination  of  the  strain  distribution  measured  at  early  low  stages  can  give  indication 
of  the  crack  initiation  location. 

5.  Crack  initiation  occurs  at  approximately  70  percent  of  the  failure  load. 

6.  Tensile  strain  levels  before  cracking  are  significantly  higher  than  those  found  under 
uniaxial  stress  conditions. 

7.  Even  though  the  Brazilian  Test  was  developed  and  has  been  useful  to  determine  the 
tensile  strength  in  concrete,  failure  does  not  start  in  the  region  of  uniform  tensile 
stress.  Thus,  the  relationship  between  tensile  strength  and  the  applied  load  may 
differ  from  that  proposed  by  Carneiro. 

8.  From  the  experimental  observations,  it  is  clear  that  stable  crack  propagation  has  an 
important  role  in  the  determination  of  the  structural  response. 

9.  Crack  initiation  location  may  be  one  of  the  reasons  why  the  post  peak  part  of  the 
load-deformation  diagram  is  difficult  to  reach.  The  effect  of  crack  opening  on  the 
diametral  displacement  is  very  small  at  the  initial  stages  of  crack  propagation. 

10.  The  reported  wedge  formation  in  the  Brazilian  Test  configuration  constitutes  a 
secondary  failure  mode  and  occurs  after  a  single  crack  has  propagated  through  the 
diameter  of  the  specimen  and  a  significant  loss  of  strength  has  taken  place.  The 
above  is  true  for  the  specimen  size  used  during  this  investigation.  A  change  of 
failure  sequence  may  occur  for  larger  specimens  (Reference  25) 


SECTION  IV 

STRAIN  FIELDS  IN  MIXED  MODE  FRACTURE  PROCESS  ZONE 


A.  INTRODUCTION 

Linear  elastic  fracture  mechanics  has  been  used  to  model  quasi-brittle  behavior  of 
cement-based  materials.  Indirect  methods  such  as  the  multicutting  technique  (Reference 
26)  and  experimental  observations  using  acoustic  emission  (Reference  27)  and 
holographic  interferometry  (Reference  28)  have  demonstrated  that  there  is  transfer  of 
tensile  stress  between  the  crack  faces.  Cohesive  crack  models  that  simulate  the  stress 
transfer  across  the  crack  and  otherwise  conventional  linear  elastic  fracture  mechanics 
have  been  used  to  simulate  Mode  I  cracking  (References  29,  7,  8  and  30). 

Since  Mode  !  tests  can  be  conducted  easily  in  laboratories  compared  to  other 
modes,  cohesive  crack  models  and  several  other  methods  have  been  proposed  to 
characterize  Mode  I  crack  propagation  in  concrete  (References  31,  32  and  33). 
However,  many  structural  applications  involve  a  mixed-mode  loading  condition.  It  is 
necessary  to  see  if  an  extension  of  the  Mode  I  model  can  be  applied  in  the  case  where 
both  tensile  and  shear  forces  associated  with  crack  opening  and  sliding  are  transmitted 
trough  the  crack  faces.  The  study  presented  in  this  section  was  intended  to  measure  the 
strain  distribution  in  the  fracture  process  zone  in  mixed  mode  specimens  in  order  to 
provide  the  information  required  to  develop  a  cohesive  crack  model  similar  to  those  for 
Mode  I. 
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B.  EXPERIMENTAL  PROGRAM 


Specimens  similar  to  the  disk  specimen  used  in  the  study  discussed  in  Section  III 
were  used.  A  0.4  inch  notch  was  cast  longitudinally  at  the  center  of  the  specimens  as 
shown  in  Figure  47.  Three  center-notched  disk  specimens  (CNDS)  were  cast  at  the  same 
time  the  unnotched  disk  specimen  for  the  Brazilian  test  were  made.  A  summary  of  the 
material  properties  is  shown  in  Table  1. 

TABLE  3  HOLOGRAPHIC  RECORDINGS.  SPECIMEN  HOLO-OO-8 
NOTCH  INCLINATION  ANGLE  B=0° 


Holographic 

Recording 

Time 

Load 

s 

lb 

N 

P  (Preload) 

350 

368 

1,656 

1 

2,500 

633 

2,849 

2 

3,000 

751 

3,380 

2* 

5,000 

670 

3,015 

2" 

7,000 

618 

2,781 

3 

8,024 

927 

4,172 

4 

8,900 

1,148 

5,166 

4* 

10,400 

1,100 

4,950 

4’* 

12,000 

1,075 

4,838 

5 

14,000 

1,457 

6,557 

6 

15,500 

1,722 

7,749 

7 

17,290 

2,105 

9,473 

3 

18,761 

2,326 

10,467 

9 

19,100 

2,355 

10,598 

10 

19,450 

2,429 

10,931 

11 

20,000 

2,400 

10,800 

12 

20,500 

2,340 

10,530 
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Figure  48  (Reference  34)  shows  how  by  rotating  the  notch  inclination  angle  (p) 
with  respect  to  the  loading  direction,  it  is  possible  to  change  the  mode  of  fracture  from 
tensile  (Mode  I  for  p=0°)  to  mixed  mode  (Tension-Shear  for  p<30°  and  Compression- 
Shear  for  p>30°). 

Both  the  mechanical  arrangement  and  the  holographic  setup  used  are  identical  to 
those  described  in  Section  III.  Diametral  displacement  was  measured  as  opposed  to  the 
crack  opening  displacement,  to  preserve  symmetry.  Preliminary  tests  indicated  that  the 
opening  force  induced  by  standard  clip  gages  is  enough  to  produce  asymmetric  cracking. 

Specimens  were  tested  at  five  different  inclination  angles  (0°,  18°,  36°,  54°  and 
72°),  but  it  was  not  possible  to  make  holographic  recordings  for  p  =  36°  during  crack 
propagation  due  to  the  brittleness  of  the  response.  The  following  two  sections  discuss 
the  results  corresponding  to  notch  inclination  angles  p  =  0°  and  p  =  54°.  In  all  cases, 
the  loading  process  was  continued  after  the  first  crack  propagated  throughout  the 
specimen,  up  to  the  limit  of  the  instrumentation  used  for  close  loop  control,  but  collapse 
of  the  specimens  was  not  reached. 

After  the  deformation  control  transducers  were  out  of  range,  the  specimens  were 
loaded  under  stroke  control,  up  to  collapse,  to  obtain  the  final  failure  mode  as  reported 
in  Figure  49. 

C.  STRAIN  FIELDS  (B  =  0°) 

Figure  50  shows  the  load  vs.  diametral  displacement  curve.  In  this  plot,  the  stars 
correspond  to  the  holographic  recordings  as  shown  in  Table  3.  The  initial  part  of  the 
curve  is  not  reliable  because  it  was  affected  by  the  vibrations  and  air  currents  induced  at 
the  start  of  the  test.  This  was  reflected  in  the  low  quality  of  the  holographic  recordings 
obtained  at  this  load  levels.  As  shown  in  Table  3,  Holograms  2  and  4  had  to  be  repeated 
until  the  desired  quality  was  obtained. 

The  shape  of  this  curve  is  similar  to  that  reported  by  Ojdrovic  and  Petroski 
(Reference  35)  for  notch-length/diameter  ratio  of  0.3.  Stable  crack  growth  occurs 
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Figure  50  Load  vs.  Diametral  Displacement  Curve.  0  =  0° 


between  Holograms  7  and  10  (2105  pounds  and  2429  pounds  respectively)  even  though 
there  is  no  macroscopic  evidence  of  nonlinearity.  After  Hologram  10,  the  first  crack 
oropagated  through  the  specimen  and  the  nonlinear  response  is  significant.  Further 
ncrease  of  the  load  after  this  point  is  possible.  A  maximum  load  of  3,040  pounds  was 
eached.  After  the  loading  sequence  was  completed,  a  second  set  of  cracks  starting 


from  the  perimeter  of  the  cylinder  have  formed,  initiating  the  wedge  formation  reported 
by  Ojdrovic. 

Figures  51  to  58  show  the  interferometric  fringes  and  corresponding  strain  fields 


Figure  51  Interferometric  Fringes.  Hologram  3 


for  four  stages  of  the  loading  process.  Figure  52  shows  the  strain  distribution  during  the 
elastic  range.  Although  there  are  variations  with  respect  to  the  linear  elastic  solution, 
there  is  indication  of  strain  concentration  in  the  region  surrounding  the  notch  tips.  Figure 
54. 

corresponds  to  crack  initiation.  At  this  load  level,  the  crack  initiated  from  the  bottom 
notch  tip.  The  strain  distribution  is  similar  around  top  and  bottom  notch  tips  because  the 
crack  is  very  short  (less  than  0.1  inch)  and  crack  initiation  at  the  top  notch  tip  is  just 
about  to  occur.  The  zones  of  strain  concentration  are  wider  and  the  strain  increase  is 
lower  compared  to  the  linear  elastic  solution.  The  linear  elastic  solution  was  obtained  by 
propagating  a  crack  in  the  finite  element  mesh  using  the  crack  trajectory  obtained  from 
the  holographic  recordings.  The  same  observation  can  be  made  during  stable  crack 
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Figure  53  Interferometric  Fringes.  Hologram  7. 


propagation  and  failure  stages  illustrated  in  Figures  56  and  58.  Also,  in  these  figures  it 
is  apparent  that  the  strain  relaxation  behind  the  crack  tip  is  complete  for  the  linear  elastic 
solution  whereas  the  experimental  strain  values  do  not  immediately  reach  zero. 

D.  STRAIN  FIELDS  (B  =  54°) 

The  load  vs.  diametral  displacement  is  shown  in  Figure  59.  A  summary  of  the 
holographic  recordings  is  in  Table  4. 

The  descending  part  of  the  load-deformation  plot  immediately  after  unstable  crack 
propagation  indicates  a  reduction  of  the  ductility  as  compared  to  the  previous  case. 
Stable  crack  growth  was  still  possible  and  occurs  between  Holograms  8  and  12  (1560 
pounds  and  2120  pounds  respectively).  A  maximum  load  of  2,307  pounds  was  reached. 
Signs  of  secondary  cracking  were  found  after  the  completion  of  the  experiment. 

Strain  distribution  and  interference  fringes  for  four  load  levels  are  shown  in  Figures 
60  to  67.  In  general,  similar  conclusions  as  those  from  the  previous  specimens  can  be 
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Figure  55  Interferometric  Fringes.  Hologram  8. 


made  regarding  the  strain  distributions  relative  to  the  crack  tip  location.  Also,  note  that 
even  at  low  load  levels,  the  strain  level  is  higher  on  the  top  half  of  the  specimen,  which 
can  be  related  to  the  mode  of  failure.  The  cracks  initiated  at  the  same  time  at  both  notch 
tips  and  initially,  the  bottom  crack  propagates  slightly.  Then,  the  top  crack  propagates 
up  to  failure. 

E.  FRACTURE  PROCESS  ZONE 

A  similar  approach  to  that  presented  in  Section  II  was  used  to  enclose  the  FPZ, 
defined  as  the  zone  of  deviation  from  linearity  exceeding  60  j Lie.  Figure  68  contains  the 
difference  in  strain  field  between  the  linear  elastic  and  the  holographic  solutions  for  notch 
inclination  angle  /3  =  0°.  Before  crack  propagation,  the  FPZ  is  larger  on  the  bottom  of 
the  specimen  which  corresponds  to  the  direction  of  maximum  crack  propagation  at  latter 
loading  stages.  The  width  of  the  fracture  process  zone  remains  approximately  constant. 
Unfortunately,  the  specimen  size  makes  unclear  whether  or  not  the  FPZ  grows  with  crack 


Figure  57  Interferometric  Fringes.  Hologram  10. 


length.  Zones  B  in  Figure  68  correspond  to  the  so-called  wake  of  the  fracture  process 
zone  (WFPZ)  and  are  attributed  to  stress  transfer  between  the  crack  faced  behind  the 
crack  tip.  Zones  B  do  not  appear  until  the  bottom  crack  has  propagated  0.2  inches  from 
the  notch  tip  which  correlates  well  with  the  model  proposed  by  Cook,  et  al.  (Reference 


29). 

If  the  same  definition  based  on  difference  in  tensile  strain  is  used,  the  FPZ  for  54°  notch 
inclination  (Figure  69)  is  smaller  than  those  obtained  for  tension  cracks  (£  =  0°  and  plate 
specimen  studied  in  Section  II)  for  short  crack  lengths.  For  longer  crack  lengths  (after 
2120  pounds  of  applied  load),  when  the  crack  stopped  curving  and  has  adopted  a 
straight  line  trajectory,  both  the  FPZ  and  WFPZ  grow  to  similar  dimensions. 

The  apparent  LEFM  stress  intensity  factors  (Figure  70)  calculated  by  Equation  (8) 

K  =  [  Kf  +  K},  11/2 
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Figure  59  Load  vs.  Diametral  Displacement  Curve.  /?  =  54° 


as  proposed  by  Jenq  (Reference  36)  follow  an  R-curve  type  of  behavior.  However,  it 
has  been  shown  in  Section  II  that  a  cohesive  crack  model  can  be  used  to  simulate  the 
stress  transfer  between  crack  faces  which  may  results  in  a  constant  (or  rather  a  narrow 
band)  of  values  for  the  stress  intensity  factor  (Reference  37).  Thus,  a  constant  fracture 
toughness  (KJ  can  be  used  to  predict  crack  propagation.  The  development  of  a 
cohesive  crack  model  for  mixed  mode  cracking  is  not  simple  because  the  complex 
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Figure  61  Holographic  Interferogram.  Hologram  8. 
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Figure  70  Apparent  Stress  Intensity  Factor^)  vs.  Crack  Extension  (a) 


1 .  Crack  initiation  occurs  at  less  than  70  percent  of  the  ultimate  load  for  both  notch 
inclination  angles.  For  notch  inclination  angle  @ =0°  the  load  can  be  increased  beyond 
the  unstable  crack  level  while  the  mixed  mode  case  is  more  brittle  (fi  =  54°)  in  the 
secondary  cracking  stage  (after  the  first  peak). 

2.  Examination  of  the  strain  distribution  at  early  loading  stages  can  give  an  indication  of 
the  crack  growth  sequence. 

3.  The  size  of  the  FPZ  defined  as  a  function  of  the  tensile  strain  distribution  can  give  an 
indication  of  the  cracking  mode.  Both  the  FPZ  and  WFPZ  are  smaller  for  mixed  mode 
cracking  than  they  are  for  mode  I. 

4.  The  mixed  mode  apparent  stress  intensity  factor  Kq  exhibits  R-curve  type  of  behavior. 


SECTION  V 


STRAIN  RATE  EFFECT 


A.  INTRODUCTION 

Experimental  results  indicate  that  the  tensile  and  compressive  strengths  of  concrete 
increase  with  increase  of  rate  of  loading  (Reference  1,  2).  Recent  investigations  have 
revealed  that  such  favorable  rate  effects  in  strength  may  be  offset  by  an  adverse  change 
in  mode  of  failure.  Reinforced  concrete  structures  designed  to  fail  in  a  ductile  flexural 
mode  at  slow  rates  of  loading  have  been  observed  to  fail  in  a  brittle  diagonal  tension- 
shear  manner  at  impact  rates  of  loading  (Reference  3). 

Knowledge  of  the  effect  of  impact  loading  on  fracture  of  concrete  is  essential  to  predict 
the  rate  dependent  failure  mechanisms  in  concrete  structures.  In  general,  the  fracture  of 
reinforced  concrete  elements  occurs  under  combined  Mode  I  and  Mode  II  loading 
conditions.  Little  information  is  available  on  the  effects  of  rates  of  loading  on  fracture 
properties  of  concrete  subjected  to  combined  tensile  and  shear  loading. 

Various  test  methods  have  been  used  to  investigate  the  effect  of  impact  loading  on 
mechanical  properties  of  concrete  in  tension,  flexure  and  compression.  Ramakrishnan, 
et.al.  (Reference  41)  used  a  conventional  drop-weight  test  to  determine  the 
performance  of  steel-fiber-reinforced  concrete.  The  test  consists  of  repeatedly  dropping 
a  known  weight  on  the  specimen.  The  impact  resistance  is  characterized  according  to 
the  number  of  strikes  required  to  induce  visible  damage.  Bhargava  and  Rehnstrom 
(42)  use  explosives  and  high-speed  photography  to  study  to  the  dynamic  cracking 


behavior  of  regular  and  polymer-modified  concrete.  The  Charpy  impact  Test  developed 
originally  for  metals  consists  of  a  pendulum-type  hammer  striking  a  specimen.  The 
fracture  energy  is  evaluated  from  the  travel  of  the  pendulum  after  impact.  However,  the 
measured  energy  also  includes  that  absorbed  by  the  testing  system  and  the  kinetic 
energy  imparted  to  the  specimen.  Both  the  drop  weight  test  and  the  Charpy  impact  Test 
have  been  modified,  adding  instrumentation  to  measure  load,  displacement  and  strain 
versus  time  responses  of  the  specimen  during  the  impact  event.  Gopalaratnam,  Shah 
and  John  (Reference  43)  developed  an  instrumented  Charpy  impact  testing  system  to 
measure  Hammer  and  support  loads  independently.  They  used  this  test  setup  to  study 
the  impact  behavior  of  steel  and  glass  fiber-reinforced  cement  composites  (Reference  44). 
For  higher  rates  of  loading,  the  split-Hopkinson  bar  test  enables  the  determination  of 
stress-strain  response  in  tension  and  compression.  Reinhardt,  et  al.,  Malvern,  et  al.  and 
Ross,  et  al.  (References  6,  5  and  4)  have  used  this  method  to  study  the  high  strain  rate 
behavior  of  cement  based  specimens  in  tension  and  compression. 

The  study  presented  in  this  section  includes  the  testing  of  the  same  center  notched 
cylinder  specimens  described  in  Section  IV  at  different  rates  of  loading.  First,  a  quasi¬ 
static  rate  of  loading  was  used  under  dosed  loop  conditions.  Reflection  photoelasticity 
and  conventional  TV  imaging  (30  frames  per  second)  was  used  to  evaluate  the  crack 
growth.  For  higher  rates  of  loading,  a  modified  Charpy  machine  and  high  speed 
photography  (10,000  frames  per  second)  were  used. 


B.  TESTING  PROCEDURE 
1.  Reflection  PtiotoelastSeity 

Many  transparent,  noncrystalline  materials  are  optically  isotropic  when  free  of  stress 
but  become  optically  anisotropic  and  display  characteristics  similar  to  crystals  when 
stressed.  This  behavior  is  known  as  temporary  double  refraction.  Materials  that  exhibit 
this  behavior  will  be  referred  to  as  photoelastic  medium  in  this  dissertation.  The  stress- 
optic  law  developed  by  Maxwell  in  1853  (Reference  45)  states  that  the  indexes  of 
refraction  are  linearly  proportional  to  the  stresses  and  followed  the  following  relationships: 
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where: 
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no 

Hi 

C1>C2 


principal  stress 
^1  +  cr2+a3 

refraction  index  in  unstressed  state 

refraction  index  in  stressed  state  associated  with  principal  direction  i 
stress-optic  coefficients 


If  a  plane  stress  condition  is  assumed,  the  stress  optic  law  from  the  previous  equation 
takes  the  form: 


Hz  -  ^1  =  (  Cr,  -  C,  )  (  a1  -  o2  )  =  C  (  Oi  -  (J2  )  (10) 


where  c  a  material  property  (usually  independent  of  the  wave  length  of  the  light)  called 
the  relative  stress-optic  coefficient. 

Now,  consider  a  photoelastic  medium  of  thickness  h  and  recall  that  the  refraction 
index  is  the  ratio  of  the  wave  velocity  in  free  space  (v0)  to  the  wave  velocity  in  the  material 
(v)  or  n=v0/v.  Thus,  the  time  (t)  required  for  passage  through  a  material  of  thickness  h 
is  given  by: 

t.*  tut 

v 

therefore,  the  distance  by  which  the  wave  trails  a  wave  traveling  on  free  space  (linear 
phase  shift)  is: 

8  -  va  t  -  h  -  ^  -  h  =  h(n  -  1  )  (12) 

v 

and  the  relative  linear  phase  shift  between  waves  traveling  along  the  principal  directions 
of  a  stressed  photoelastic  medium  is: 
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which  corresponds  to  the  relative  angular  phase  shift  given  by: 


A  =  —  8 
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(14) 


In  engineering  practice,  this  expression  is  usually  expressed  in  the  following  form: 
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i 


(15) 


or  in  terms  of  principal  strains  (if  the  photoelastic  medium  is  perfectly  elastic): 


s 


i 


(16) 


where  N  is  the  relative  phase  shift  in  terms  of  complete  cycles  of  retardation  and  fa  and 
fs  are  known  as  the  material  fringe  values  in  terms  of  stress  and  strain  respectively  and 
can  be  established  by  calibration. 

A  polariscope  is  used  in  order  to  measure  the  phase  shift  (N)  at  each  point  on  a 
stressed  model.  A  dark  field  circular  polariscope  (Figure  71)  was  used  during  this 
investigation.  A  detailed  explanation  of  the  operation  of  this  and  other  types  of 
polariscope  was  shown  in  Reference  45. 

The  objective  of  the  polariscope  is  to  transform  relative  phase  shifts  into  observable 
intensity  differences.  If  the  stressed  photoelastic  medium  shown  in  Figure  71  is 
illuminated  through  the  polariscope  by  a  light  beam  of  intensity  k  and  frequency  u,  the 
light  wave  emerging  from  the  polariscope  is: 
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Figure  71  Circular  Polariscope  (Dark  Field) 

E„  =  /rsin^  sin|«f  +2a  -  -|-j  (17) 

Since  the  intensity  of  light  is  only  proportional  to  the  square  of  the  amplitude  of  the 
light  wave,  the  intensity  of  the  light  beam  emerging  from  the  circular  polariscope  is  a 
function  only  of  the  angular  phase  shift  (independent  of  the  frequency,  o>  and  the  phase, 
a)  and  is  given  by: 


99 


When  A/2  =n7r,  where  n=0,1,2,3,...,  the  intensity  value  is  zero  (extinction).  A  fringe  is 
formed  of  the  loci  of  all  the  point  in  the  model  with  the  same  order  of  extinction  (n).  This 
fringe  pattern  is  called  the  isochromatic  fringe  pattern.  The  order  of  extinction  n,  depends 
on  both  the  principal  stress  difference  and  the  wave  length  used  to  illuminate  the  model 
(Equation  V.7).  When  the  model  is  illuminated  with  white  light  (linear  combination  of  light 
waves  of  different  wave  length),  the  isochromatic  fringe  pattern  appears  as  a  colored 
fringe  pattern  and  there  is  a  relationship  between  fringe  color  and  relative  phase  shift. 
The  polariscope  described  above  is  called  a  dark  field  polariscope  because  the 
unstressed  material  appears  dark.  If  the  analyzer  plate  is  rotated  SO  degrees  such  that 
the  axis  of  polarization  is  parallel  to  that  of  the  polarizer  plate  (light  field  polariscope),  the 
unstressed  material  would  appear  white.  In  this  investigation  a  dark  field  arrangement 
was  selected  for  quasi-static  loading  and  a  light  field  polariscope  for  impact  loading 
because  they  were  best  suited  for  the  corresponding  imaging  method. 

In  order  to  study  opaque  objects,  the  method  of  birefringent  coatings  or  reflection 
photoelasticity  is  used.  The  photoelastic  medium  is  a  thin  coating  which  is  bonded  to  the 
surface  of  the  specimen.  Figure  72  shows  the  reflection  polariscope  used  for  this  study. 
The  operation  of  the  polariscope  is  very  similar  to  that  of  a  conventional  circular 
polariscope.  The  bonding  agent  between  the  photoelastic  medium  and  the  specimen  is 
highly  reflective  and  provides  a  mirror  like  surface.  Since  the  light  travels  through  the 
coating  in  both  the  incident  and  the  reflected  direction,  a  modified  thickness  h’=2h  should 


be  used  for  the  analysis. 


Specimen 


Figure  72  Circular  Reflection  Polariscope 


of  the  specimen  to  preserve  symmetry.  Specimens  were  tested  using  the  same 
mechanical  setup  and  instrumentation  described  in  Section  III  up  to  the  limit  of  the 
instrumentation  used  for  closed  loop  control.  A  Panasonic  S-VHS  TV  camera  was  used 
to  record  the  isochromatic  fringes  continuously. 

3.  Instrumented  Impact  Loading 

The  instrumented  Charpy  machine  used  by  Gopalaratnam,  Shah  and  John  (Reference 
43)  (Figure  73)  was  modified  to  accommodate  CNDSs.  A  modified  instrumented  striker 
(tup),  instrumented  anvil,  a  specimen  support  and  an  anchor  were  fabricated.  Figures  74 
to  76  show  the  instrumented  striker,  instrumented  tup  and  general  assembly  for  he 
modified  Charpy  machine.  A  full  Wheatstone  bridge  was  connected  on  both  tup  and  anvil 
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Figure  73  Instrumented  Charpy  Machine 


sensitive  only  to  axial  load.  A  Tektronix  2232  model  oscilloscope  was  used  for  data 
acquisition  at  a  rate  of  250  KHz.  Both  the  striker  and  the  anvil  were  calibrated  at  static 
rate  of  loading  in  a  20  kip  capacity  MTS  load  frame  under  load  control  up  to  10,000  lb. 
Dynamic  calibration  was  achieved  by  placing  a  dynamic  load  cell  in  place  of  the 
specimen.  Typical  calibration  results  are  shown  in  Figure  77.  In  order  to  estimate  the 
inertial  effect  dummy  specimens  were  tested  and  the  responses  from  the  tup  and  anvil 
were  compared  as  shown  in  Figure  78.  It  was  found  that  both  peak  load  and  arrival  time 
differences  are  small  (5  percent)  compared  to  the  overall  response  and  inertia  effects  can 
be  neglected. 
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Figure  74  Instrumented  Striker  (Tup) 
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Figure  76  Modified  Charpy  Machine  Striker-Support  Assembly 


For  evaluation  of  crack  propagation,  the  same  photoelastic  coating  described  in 
Section  V. 2.2  was  glued  to  front  and  back  faces  of  the  specimens.  A  lens  was 
incorporated  in  the  polariscope  arrangement  (Figure  79)  to  increase  the  light  intensity  to 
the  level  required  for  high-speed  photography.  Halogen  light  was  used  for  illumination. 
A  high-speed  camera  Model  100  (Photographic  Analysis  Company)  was  used  for 
recording  of  the  isochromatic  fringes  (Figure  80).  450-foot  rolls  of  SO-259  Kodak  film 
(400  ASA)  for  high  speed  photography  were  used.  The  high  speed  camera  was  used  to 
trigger  the  overall  system.  After  a  time  delay  required  to  bring  the  film  up  to  speed,  the 
camera  closes  an  electric  circuit  which  activates  a  relay  to  release  to  pendulum  of  the 
Charpy  machine  (Figure  81).  Data  acquisition  was  triggered  by  the  impact  between  the 
striker  and  the  specimen.  A  second  2232  Model  Tektronix  oscilloscope  was  used  at  low 
frequency  (5  KHz)  to  measure  the  time  difference  between  the  release  of  the  pendulum 
and  impact  with  the  specimen. 
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Figure  78  Dynamic  Calibration.  (Inertia  Effect) 

In  the  load-deformation  curve  for  the  unnotched  specimen  shown  in  Figure  83,  the 
stars  indicate  the  location  of  photographic  recordings.  As  discussed  earlier,  continuous 
imaging  at  30  images  per  second  were  recorded  using  a  S-VHS  tape  recorder.  The  most 
significant  stages  of  crack  propagation  are  shown  in  Figures  84  and  85.  Note  that  the 
failure  sequence  follows  closely  the  holographic  observations  presented  in  Section  III. 5. 
The  images  in  Figure  84  and  85  lead  to  the  conclussion  that  the  more  sensitive  (0.125 
inch  thickness)  photoelastic  coating  used  in  this  experiment  is  sensitive  to  deformations 
along  the  crack  trajectory,  but  measurement  of  the  far-field  deformations  is  not  possible. 
To  reduce  the  ’reinforcing’  effect  of  the  photoelastic  coating,  a  thinner,  less  sensitive 
coating  was  used  for  all  the  subsequent  experiments. 
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Figure  80  High  Speed  Photography  Arrangement  for  Instrumented  Impact 


Figure  86  shows  the  load  deformation  curve  for  a  specimen  with  notch  inclination 
angle  f3  =  0°.  Note  that  at  Point  B,  there  is  a  first  peak  characteristic  of  crack  initiation 
and  propagation  based  on  the  holographic  results  reported  earlier.  However,  the 
isochromatic  fringe  patterns  shown  in  Figures  87  and  88  indicate  that  cracking  starts  at 
Point  C.  Between  Points  C  and  G  the  main  crack  grows  through  the  specimen  and  a 
second  peak  is  observed  in  the  load-deformation  curve.  The  above  phenomenon  can  be 
attributed  to  asymmetric  failure.  Previous  calibration  experiments  had  shown  that  it  is 
possible  for  the  crack  to  initiate  and  propagate  on  one  face  of  the  specimen  before 
cracking  starts  on  the  opposite  face  due  to  slight  eccentricity  of  the  load  or  specimen 
imperfections.  A  difference  in  behavior  between  the  holographic  test  reported  in  Section 
II  and  this  experiment  is  the  significant  additional  load  capacity  after  crack  propagation. 
After  cracking,  a  secondary  failure  mechanism  develops  which  results  in  the  failure  wedge 
reported  in  Section  IV.  However,  even  after  large  crack  openings,  the  photoelastic 
coating  provides  a  relatively  strong  joint  between  the  crack  faces  as  compared  to  the 
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Figure  81  Pendulum  Release  Mechanism 


uncoated  specimens.  Final  failure  occurs  at  2.5x10‘3  inch  of  diametral  displacement, 
which  corresponds  to  the  wedge  formation. 

Figures  89  to  91  show  the  load-deformation  curve  and  isochromatic  fringe  patterns  for 
a  specimen  with  notch  inclination  angle  fS  =  54°.  Symmetric  failure  was  observed  in  this 
specimen  and  both  the  load-deformation  curve  and  the  crack  growth  agree  with  the 
observation  obtained  using  laser  holographic  interferometry.  The  first  crack  propagates 
through  the  specimen  between  Points  B  and  G  located  before  the  first  peak  in  the  load- 
deformation  curve.  The  second  peak  reaches  a  slightly  higher  load  value  than  the  first 
peak  (as  opposed  to  the  lower  load  level  in  the  holographic  experiment).  The  top  of  the 
second  peak  can  be  correlated  to  the  initiation  of  the  secondary  crack  at  Point  K. 
Between  the  two  peaks  the  physical  mechanism  is  characterized  by  opening  of  the  first 
crack  and  the  reinforcing  effect  of  the  photoelastic  coating  becomes  important. 

After  the  secondary  cracks  have  grown,  there  is  a  plateau  in  the  load-deformation  curve 
at  collapse  of  the  specimen  occurs  at  point  N  after  the  photoelastic  coating  has 
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Figure  82  Load  Deformation  Curves  for  Specimens  Tested  at  Quasi-Static  Loading. 


debonded  from  the  specimen. 

Post-failure  examination  of  the  specimen  indicates  that  the  concentration  of 
isochromatic  fringes  follows  closely  the  crack  trajectory  on  the  surface  of  the  specimen. 
Thus,  the  possibility  of  local  debonding  of  the  photoelastic  coating  before  cracking  should 
be  eliminated. 

In  general,  an  acceptable  agreement  was  found  between  holographic  and 
photoelasticity  observations  regarding  crack  formation  and  propagation  at  slow  rates  of 
loading.  The  modes  of  failure  were  not  modified  by  the  existence  of  the  photoelastic 
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Figure  83  Load  Deformation  Curve.  Unnotched  Specimen.  Quasi-Static  Loading. 

coating  and,  although  some  reinforcing  effects  were  observed  the  structural  behavior  of 
the  specimen  showed  the  same  trends  for  different  notch  inclination  angles. 

2,  Impact  Loading  (strain  rate  1  to  10  per  sec) 

A  comparison  of  the  load  vs.  recordings  at  impact  loading  for  various  notch  inclination 
angles  is  shown  in  Figure  92.  In  general,  the  trends  agree  with  results  obtained  at  low 
rates  of  loading.  For  inclination  angles  /3=0°  and  /3  =  18°,  there  is  evidence  of  ductile 
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Figure  84  Isochromatic  Fringes.  Unnotched  Specimen.  Quasi-Static  Loading. 
(Expossures  A,B  and  C) 
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Figure  85  Isochromatic  Fringes.  Unnotched  Specimen.  Quasi-Static  Loading. 
(Expossures  D,E  and  F) 
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Figure  86  Load-Deformation  Curve.  0=0°.  Quasi-Static  Loading. 


behavior  represented  as  a  short  plateau  at  the  peak  load  while  for  higher  inclination 
angles  the  response  is  quite  brittle.  The  unnotched  specimen  shows  a  sharp  unloading 
after  the  peak  load  followed  by  a  subsequent  region  of  constant  load  over  a  relatively  long 
period  of  time.  Two  unnotched  specimen  were  tested  under  impact  loading  using  high 
speed  photography.  Unfortunately,  none  of  the  recordings  was  successful  in  observing 
crack  growth  for  the  unnotched  specimen. 

The  impact  energy  was  selected  such  that  the  failure  mode  remained  the  same  as  that 
observed  during  the  slow  rate  experiments.  Therefore,  the  observed  macroscopic  rate 
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Figure  88  Isochromatic  Fringes,  p  =  0°.  Quasi-Static  Loading. 
(Expossures  G  and  L) 
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Figure  89  Load-Deformation  Curve.  p  -  54°.  Quasi-Static  Loading. 

effects  can  be  related  only  to  the  rate  effect  on  crack  initiation  and  crack  growth.  Higher 
rates  of  loading  usually  result  in  fragmentation  of  the  specimen  and  different  failure 
mechanisms. 

Figure  93  shows  the  load  record  for  a  specimen  with  notch  inclination  angle  (3=0°. 
An  isochromatic  fringe  pattern  was  recorded  at  constant  intervals  of  0.1  ms.  however, 
the  timing  of  the  photographic  exposures  can  only  be  determined  within  0.05  ms  (i.e. 
photographic  exposures  could  be  shifted  forward  or  backward  by  0.05  ms  in  Figure  93 
but  the  interval  between  exposures  remains  0.1  ms).  Figures  94  and  95  show  the 
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Figure  92  Load  History  Record  for  Various  Notch  Inclination  Angles  at  Impact 
Loading. 

sequence  of  isochromatic  fringe  pattern  recordings  corresponding  for  the  same 
specimen.  Unfortunately,  the  pictures  shown  in  this  dissertation  are  not  as  clear  as  the 
original  film  due  to  the  reproduction  process  but  some  observation  can  still  be  made.  The 
failure  sequence  agrees  with  both  the  holographic  and  photoelasticity  observations  at  low 
rates  of  loading.  Point  C,  shows  a  concentration  of  deformation  along  the  crack  path  but 
it  is  unclear  whether  crack  has  initiated  or  not.  At  Point  D  there  is  a  crack  on  each  tip 
of  the  notch  and  started  to  propagate  towards  the  points  of  application  of  the  load  and 
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appears  to  be  an  unloading  region  in  the  load  vs.  time  record.  This  corresponds  to  first 

peak  is  found  at  po 

secondary  crack  initiates  (Figure  95).  Postfailure  examination  of  the  specimens  showed 
that  the  crack  patterns  correspond  to  the  isochromatic  fringe  patterns.  Debonding  of  the 
photoelastic  coating  was  observed.  Further  examination  of  the  failed  specimens  showed 


Figure  94  Isochromatic  Fringe  Patterns.  /3  =  0°.  Impact  Loading. 
(Expossures  C  and  D) 
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Figure  95  isochromatic  Fringe  Patterns.  /3  =  0°.  Impact  Loading. 
(Expossures  E  and  G) 
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back  face  of  the  photoelastic  medium. 

The  load  history  for  a  specimen  with  notch  inclination  angle  p  =  54°  is  shown  in 
Figure  96.  Figures  97  and  98  show  the  sequence  of  isochromatic  fringe  pattern 


Figure  96  Load  Record  and  Photographic  Exposures,  p  =  54°.  Impact  Loading. 

recordings  corresponding  to  this  specimen.  The  isochromatic  fringe  pattern  at  point  B 
shows  a  concentration  of  strain  around  both  notch  tips  and  the  first  crack  are  detected 
at  point  C.  At  point  D  the  crack  has  propagated  to  the  perimeter  of  the  specimen  and 
final  failure  is  accompanied  by  formation  of  secondary  cracks  and  debonding  of  the 
photoelastic  coating  at  point  F. 
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Figure  97  Isochromatic  Fringe  Patterns,  p  =  54°.  impact  Loading. 
(Expossures  B  and  C) 
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Figure  98  Isochromatic  Fringe  Patterns.  0  =  54°.  Impact  Loading. 
(Expossures  D  and  F) 


127 


retardation  of  the  crack  growth  compared  to  the  loading  sequence.  Thus,  it  can  b 
concluded  that  for  the  rates  of  loading  used  during  this  investigation  failure  i 
characterized  by  stable  crack  growth.  Crack  propagation  velocity  at  this  load  rat 


figure  it  can  be  concluded  that  the  increase  in  apparent  fracture  toughness  for  a  tension 
crack  {p =0°)  is  higher  than  that  for  a  shear  crack  (p=36°).  Also,  note  that  between  notch 
inclination  angles  >3  =  18°  and  36°  the  relative  Kq  is  minimum,  which  corresponds  to  the 
maximum  shear  component.  In  Figure  99  the  values  corresponding  to  p  =  72°  should 
not  be  considered  very  reliable  because  the  cracks  did  not  start  from  the  notch  tips  as 
expected.  Also,  note  in  Figure  100  that  for  the  rates  of  loading  used  during  this 
investigation  the  failure  modes  remain  essentially  the  same.  Thus,  any  variation  in  the 
macroscopic  response  should  be  related  to  change  in  fracture  toughness. 

D.  CONCLUSIONS 

1 .  Reflection  photoelasticity  can  be  used  to  determine  crack  initiation  and  crack  growth 
at  both  quasi-static  and  impact  rates  of  loading.  However,  the  thicker,  more  sensitive 
coatings  are  not  sensitive  enough  to  provide  information  regarding  the  state  of 
deformation  away  from  the  crack  path. 

2.  The  reinforcing  effect  of  the  photoelastic  coating  is  evident  in  the  post-peak  region  but 
it  is  not  large  enough  to  modify  the  failure  mode. 

3.  There  is  good  agreement  between  the  observations  made  using  holographic 
interferometry  and  reflection  photoelasticity  regarding  crack  initiation,  crack  growth  and 
failure  mechanisms  at  quasi-static  loading. 

4.  At  both  quasi-static  and  impact  loading,  the  state  of  cracking  can  be  correlated  to  the 
loading  stage  of  the  specimen.  For  all  notch  inclination  angles,  there  are  three  loading 
regions  limited  by  crack  growth,  crack  opening,  and  secondary  cracking.  At  the  rates 
of  loading  used  during  this  investigation,  failure  appears  to  be  characterized  by  stable 
crack  growth. 

5.  The  increase  in  fracture  toughness  at  high  rates  of  loading  is  affected  by  the  mode  of 
fracture.  The  increase  of  fracture  toughness  for  tensile  cracking  is  approximately 
double  the  corresponding  to  mode  II  cracking. 
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Figure  100  Failure  Modes  at  Quasi-Static  and  Impact  Loading 


SECTION  VI 


CONCLUSIONS 


A.  HOLOGRAPHIC  INTERFEROMETRY 

Sandwich  hologram  interferometry  with  multiple  sensitivity  vectors  can  be  used  to 
measure  crack  opening  displacements  and  displacement  fields  in  two  and  three 
dimensions  under  closed  loop  testing  conditions.  When  holographic  interferometry  is 
used  to  measure  Mode  I  crack  opening  displacements  the  effect  of  out  of  plane  motion 
should  be  taken  into  account. 

Digital  image  analysis  facility  provided  a  faster,  more  consistent  and  more  accurate 
method  of  fringe  count  than  the  manual  method.  Image  processing  was  essential  for  the 
elimination  of  rigid  body  motion  without  loss  of  fringe  order  accuracy. 

B.  REFLECTION  PHOTOELASTICITY 

Reflection  photoelasticity  can  be  used  to  determine  crack  initiation  and  crack  growth 
at  both  quasi-static  and  impact  rates  of  loading.  Even  the  thicker,  more  sensitive  coatings 
are  not  sensitive  enough  to  provide  information  regarding  the  state  of  deformation  away 
from  the  crack  path.  There  is  good  agreement  between  the  observations  made  using 
holographic  interferometry  and  reflection  photoelasticity  regarding  crack  initiation,  crack 
growth  and  failure  mechanisms  at  quasi-static  loading. 
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C.  MODE  i  FRACTURE  PROCESS  ZOWE 

There  is  a  region  of  high  tensile  strain  that  moves  with  the  tip  of  the  propagating 
crack.  The  material  behind  the  crack  tip  experiences  strain  relaxation.  However,  the 
gradual  nature  of  the  strain  relaxation  demonstrates  the  existence  of  tensile  force 
transmitted  through  the  crack  faces. 

Zones  of  nonlinear  behavior  (FPZ)  can  be  located  by  computing  strain  field  deviations 
form  the  linear  elastic  solution.  A  bilinear  closing  pressure  vs.  COD  cohesive  crack  model 
based  on  crack  profile  measurements  has  been  proposed  to  predict  both  crack  profiles 
and  strain  field  in  mode  I  cracking,  fully  characterizing  the  material  behavior. 


D.  BRAZSLIAM  TEST 

Examination  of  the  strain  distribution  measured  at  early  low  stages  can  give  indication 
of. the  crack  initiation  location.  After  crack  initiation,  the  distribution  of  measured  strains 
can  give  an  indication  of  the  crack  growth  sequence. 

Although  the  Brazilian  Test  was  developed  and  has  been  useful  to  determine  the 
tensile  strength  in  concrete,  failure  does  not  start  in  the  region  of  uniform  tensile  stress. 
Tedesco,  et  al  (Reference  24)  made  a  similar  observation  under  high  rate  of  loading. 
Thus,  the  relationship  between  tensile  strength  and  the  applied  load  may  differ  from  that 
proposed  by  Carneiro.  From  the  experimental  observations,  it  is  clear  that  stable  crack 
propagation  has  an  important  role  in  the  determination  of  the  structural  response. 

The  reported  wedge  formation  in  the  Brazilian  Test  configuration  constitutes  a 
secondary  failure  mode  and  occurs  after  a  single  crack  has  propagated  through  the 
diameter  of  the  specimen  and  a  significant  loss  of  strength  has  taken  place. 


E.  MIXED  MODE  FRACTURE  PROCESS  ZO^SE 

The  level  of  observed  strains  at  crack  initiation  is  considerably  higher  for  the  diametral 
compression  disk  specimen  (biaxial  state  of  stress)  as  compared  to  the  compact  tension 
specimen  (uniaxial  state  of  stress). 
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The  size  of  the  FPZ  defined  as  a  function  of  the  tensile  strain  distribution  can  give  an 
indication  of  the  cracking  mode.  Both  the  FPZ  and  WFPZ  are  smaller  for  mixed  mode 
cracking  than  they  are  for  Mode  I. 

The  mixed  mode  apparent  stress  intensity  factor  Kq  exhibits  R-curve  type  of  behavior. 


F.  RATE  EFFECT 

The  observed  mode  of  failure  at  both  quasi-static  and  impact  loading  remained 
essentially  the  same  which  indicates  that  the  variation  in  response  due  to  change  in  the 
rate  of  loading  depends  on  crack  formation  and  growth. 

For  the  specimen  used  during  this  investigation,  after  initial  crack  growth,  the  crack 
continued  to  propagate  in  Mode  I  regardless  of  the  crack  initiation  mode.  Collapse  of  the 
specimen  was  characterized  by  a  secondary  crack  formation.  The  ultimate  load  for 
various  inclination  angles  did  not  dependent  on  the  mode  of  crack  initiation  alone.  Thus, 
comparison  of  the  ultimate  load  does  not  reflect  mode  of  failure  dependence.  The  use 
of  high-speed  photography  allowed  evaluation  of  the  fracture  toughness  at  crack  initiation. 

The  increase  in  apparent  fracture  toughness  at  crack  initiation  for  high  rates  of  loading 
is  affected  by  the  mode  of  failure.  The  increase  of  apparent  fracture  toughness  for  tensile 
cracking  is  approximately  double  the  corresponding  value  for  mode  II  cracking. 
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